Neutrino signal from 7-ray loud binaries powered by high energy protons 
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We present a hadronic model of activity for Galactic 7-ray-loud binaries, in which the multi- 
TeV neutrino flux from the source can be much higher and/or harder than the detected TeV 7-ray 
flux. This is related to the fact that most neutrinos are produced in pp interactions close to 
the bright massive star, in a region optically thick for the TeV 7-rays. Considering the speciflc 
example of LS I +61° 303, we derive upper bounds for neutrino fluxes from various proton injection 
spectra compatible with the observed multi-wavelength spectrum. At this upper level of neutrino 
emission, we demonstrate that ICECUBE will not only detect this source at 5o" C.L. after one 
year of operation, but, after 3 years of exposure, will also collect a sample marginally sufficient to 
constrain the spectral characteristics of the neutrino signal, directly related to the underlying source 
acceleration mechanisms. 



I. INTRODUCTION 

The recently discovered "7-ray-loud" binaries (GRLB) 
form a new sub-class of Galactic binary star systems 
which emit GeV-TeV 7-rays P, i, i, i, H. These are 
high mass X-ray binaries (HMXRB) composed of a com- 
pact object (a black hole or a neutron star) orbiting a 
massive star. The detection of 7-rays with energies up to 
10 TeV from these systems shows that certain HMXRBs 
host powerful particle accelerators producing electrons 
and/or protons with energies above 10 TeV. 

At the moment it is not clear which physical process 
leads to the very high energy (VHE) particle acceleration 
and whether the VHE particle acceleration is a generic 
feature of the HMXRBs or the result of specific physical 
conditions in a restricted HMXRB sub-class. It is pos- 
sible that particle acceleration is taking place in a large 
number of X-ray binaries, but the 7-rays can be detected 
only from the objects with preferred orientations w.r.t. 
the Hue of sight (e.g. with jets pointing toward the ob- 
server) . Alternatively, it may be that the conventional 
accretion-powered binaries are not capable of accelerat- 
ing particles and that the 7-ray activity of a binary is 
powered by a different mechanism (see e.g. 0, 0)- 

The theoretical models of 7-ray activity of HMXRBs 
0, i, i, [3, m, E Q [11 all assume 7-ray emission 
from the interaction of a relativistic outflow from the 
compact object (jet from a black hole, or wide angle wind 
from a pulsar) with the wind and radiation emitted by 
the companion massive star. The basic properties of the 
outflows, such as the composition (e"'"e^ pairs or electron- 
nuclei plasma), anisotropy (a coUimated jet or a wide 
angle outflow) etc. are as yet poorly constrained by the 
data. 

In one of the three persistent (periodic) 7-ray-loud bi- 
naries, PSR B1259-63, the relativistic outflow is known 
to be produced by a young pulsar. In principle, a simi- 
lar mechanism could power the activity of other sources 
(except for Cyg X-1), although direct proof of the pres- 
ence of the young pulsar in these systems is not possible: 
the radio emission from the compact object is free-free 



absorbed. Recent puzzling detection of a short soft 7-ray 
flare from one of the 7-ray-loud binaries, LSI 4-61 303 
p^ . may indicate the presence of a neutron star in the 
system. The outflows from LSI +61 303 and LS 5039, 
extending to the distances ~ 10^^ cm (far beyond the bi- 
nary orbit) are revealed by the radio observations (l7l.lT§|. 
The nature of the outflows (a coUimated jet or a wide- 
angle outflow) is still debated. 

With the exception of Cyg X-1, all the known GRLBs 
have similar spectral energy distributions (SED), peaking 
in the MeV-GeV energy band. The physical mechanism 
producing the MeV-GeV bump in the spectra is, however, 
not clear. The 7-ray emission from GRLBs is supposed 
to come from internal and/or external shocks formed in 
the relativistic outflow either as a result of the devel- 
opment of intrinsic instabilities or through interactions 
with the stellar wind of the massive star. The MeV-GeV 
emission can be the synchrotron emission from electrons 
with energies much above TeV @, produced locally 
at the shock, or else, be produced via inverse Compton 
scattering of the UV thermal emission from the massive 
star by electrons of the energies 10 McV 0,[ii,iil. 

The available multi-wavelength data do not constrain 
the composition of the relativistic wind from the com- 
pact object. On the one hand, the multi-TeV or 10 MeV 
electrons, responsible for the production of the MeV-GeV 
bump in the SED, could be injected into the shock re- 
gion from the e"'"e~ pairs loaded wind. On the other, 
these electrons could be secondary particles produced in 
e.g. proton-proton collisions, if the relativistic wind is 
proton-loaded. The only direct way to test if relativistic 
protons are present in the 7-ray emission region would 
be the detection of multi-TeV neutrinos. 

The possibility of detection of neutrinos from GRLBs 
by a km'^-class neutrino telescopes was considered in sev- 
eral references [H, [13, El 12, J3i, [H] . If the sources are 
assumed to be transparent to the TeV 7-rays, the esti- 
mation of the neutrino flux is straightworward, given a 
known source 7-ray flux and spectrum. The assumption 
of source transparency was adopted e.g. in the estimates 
of the number of detected neutrinos as in Ref. [231 . How- 
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ever, the TeV 7-ray flux from the GRLBs can be signifi- 
cantly attenuated by the pair production on the UV pho- 
ton field of the massive star [13, [Hi. In addition, if the 
7-rays and neutrinos are produced close to the compact 
object, the 7-ray flux can be further suppressed by pair 
production on the soft photons emitted by the accretion 
flow ^d^. The derivation of the estimate of the neutrino 
flux and spectral characteristics based on the observed 
TcV 7-ray emission is inconclusive due to the uncertain 
attenuation of the 7-ray flux in the TeV band (see [l^l 
for s specific discussion of LS I -1-61° 303). 

In the absence of a direct relation between the charac- 
teristics of the observed TeV 7-ray and the neutrino emis- 
sion from a GRLB, the only way to constrain possible 
neutrino signals from the source is via the detailed mod- 
elling of the broad-band spectrum of the source within 
the hadronic model of activity. The idea is that the pp 
interactions, which result in the production of neutrinos, 
also result in the production of e~''e~ pairs and the sub- 
sequent release of their energy via synchrotron, inverse 
Compton and Bremsstrahlung emission. The electromag- 
netic emission from the secondary e+c^ pairs is readily 
detectable. The total power released in the pp interac- 
tions determines the overall luminosity of the emission 
from the secondary pairs. The known electromagnetic 
luminosity and broad-band spectral characteristics of the 
source can be used to constrain the power released in pp 
interactions as well as the spectrum of the primary high 
energy protons. 

In the following we develop the hadronic model of ac- 
tivity of GRLBs and derive the constraints on the spec- 
trum and overall luminosity of neutrino emission from 
the analysis of the broad-band spectral characteristics 
of GRLBs. Although the following discussion is generi- 
cally applicable for the GRLBs as a class, we concentrate 
on the particular example of the LS I -1-61° 303 system, 
because it is the only known persistent GRLB in the 
Northern hemisphere, available for observations with the 
ICECUBE neutrino telescope [1^ . The existing observa- 
tional data are found to be consistent with a possible very 
strong neutrino emission from LS I -1-61° 303 (with a flux 
at the level of 10~^° erg/cm^s), close to the best reported 
AMANDA upper limit, «>(1.6 TeV < < 2.5 PeV) < 
1.26 X 10-1° cm-i s-i TeV-i) ^ (a more recent ref- 
erence [13 combining more data reported a slightly de- 
graded upper limit for LS I -1-61° 303 over approximately 
the same range of energy, while the sensitivity had been 
improved by roughly a factor two; this however translates 
a non significant but large background fluctuation at the 
LS I -1-61° 303 source location). Moreover, the assump- 
tion of an almost arbitrarily hard neutrino spectrum (e.g. 
a powerlaw dN,j/dE ^ E'^" with index T^, ~ 0) is not 
ruled out by the available multi- wavelength observational 
data. This means that, in principle, the soure could be 
readily detectable with ICECUBE. 

The plan of the paper is as follows: In section [Til wc 
describe the general features of the hadronic model of ac- 
tivity of GRLBs. In particular, we stress that when the 



massive companion star is a Be type star, the presence 
of a dense equatorial decretion disk around the massive 
star can boost the pp interaction rate. At the same time, 
the 7-ray emission from the pp interactions in the disk 
would be strongly suppressed, because of the large den- 
sity of the soft photon background in the direct vicinity 
of the star. In section IIIII wc perform the detailed nu- 
merical modeling of the broad-band spectra of GRLBs in 
the hadronic model and show, on the particular exam- 
ple of LS I -1-61° 303, that the model provides a suitable 
explanation of the typical shape of the GRLB SEDs. Sat- 
isfactory fits to the observed SEDs can be achieved un- 
der quite arbitrary assumptions about the shape of the 
initial high energy proton spectrum. In particular, the 
initial spectra as different as an E~'^ power law with a 
high energy cut-off and a monochromatic proton spec- 
trum peaked at '-^PeV energy could both explain the ob- 
served X-ray-to-7-ray spectrum of the source. Finally, in 
section IIVI we work out the predictions for the detection 
of neutrinos from LS I -^61° 303 with ICECUBE. We 
show that in an optimistic scenario, when the anisotropy 
of the neutrino emission of the source does not result 
in a suppression of the neutrino flux in the direction of 
the Earth, LS I -|-61° 303 should be detectable within a 
year of exposure. Nevertheless, the spectral character- 
istics of the emission can only be marginally delineated 
after three years of exposure, given the wash out of the 
original neutrino spectrum from the measurement of the 
muon spectrum. 



II. HADRONIC MODEL OF 7-RAY ACTIVITY 
A. Origin of high energy protons 

In the hadronic model, the primary source of the sys- 
tem's high energy activity are high energy protons. The 
presence of the high energy protons in relativistic out- 
flows from compact objects (stellar mass and supermas- 
sive black holes, neutron stars) is usually difficult to de- 
tect, because of their very low energy loss rates. For ex- 
ample, in the case of relativistic winds ejected by young 
pulsars, it is possible that most of the wind power is car- 
ried by the high energy protons or ions, but the presence 
of protons/ions in the wind can be established only in- 
directly, because most of the radiation detected from a 
nebula powered by the pulsar wind is emitted by elec- 
trons accelerated at the shock whose properties are de- 
termined by the parameters of the proton/ion component 
of the pulsar wind [251] . 7-ray emission initiated by in- 
teractions of the pulsar wind protons with the ambient 
medium can be detected only if the density of the medium 
is high enough ^2^, . 

GRLBs provide a unique opportunity to "trace" the 
presence of protons/ions in the relativistic outflows gen- 
erated by compact objects. The dense matter and radia- 
tion environment created by the companion massive star 
provides abundant target material for the protons in the 
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relativistic outflow. Interactions of high energy protons 
with the ambient matter and radiation fields, created by 
the presence of a bright massive companion star, lead 
to the production and subsequent decays of pions. This 
results in the emission of neutrinos and 7-rays from the 
source and to the deposition of e^e~ pairs throughout the 
proton interaction region. Radiative cooling of the sec- 
ondary e"'"e~ pairs leads to the broad-band synchrotron 
and inverse Compton emission from the source. 

If the primary high energy protons are produced close 
to the compact object, the maximal energies of particles 
accelerated in the region of the size in magnetic field 
Be, can be estimated as 



Ep ~ KeBcRc ^ lO^^K 



Br 



10^ G 



Rc 



106 



eV (1) 



The parameter k < 1 characterizes the acceleration pro- 
cess efficiency. For example, in the case when the com- 
pact object is a rotation-powered neutron star, B ^ 
10^2 G, while K ~ ^iRl/c^ 4 X 10"^ [Pc/0.1 s]"^ 
with i?c,^^c,^c being respectively the neutron star ra- 
dius, angular velocity and rotation period. Depending on 
the particular acceleration mechanism, the proton energy 
spectrum below this maximal energy can range from the 
conventional power law shape with spectral index close to 
Fp ~ 2 (acceleration in non-relativistic shocks), to almost 
monochromatic spectra (if the acceleration proceeds in 
large scale electric fields in the magnetic reconnection re- 
gions of accretion flow, or in the vacuum gaps in black 
hole or pulsar magnetosphere). 

Otherwise, protons can be accelerated via shock ac- 
celeration in an extended region of the size i?oxt with 
magnetic field i^oxt in a jet or in a shocked relativistic 
wind emitted by the compact object. In this case one ex- 
pects a conventional E~'^ type spectrum with a cut-off at 
the energy at which the Larmor radius of the high energy 
particles becomes comparable to the size of the system 
Ep ~ 0.3 X 1015 [Bext/1 G] [i?cxt/10i2 cm] cV. 



B. Interactions of high energy protons 

Efficient interaction of high energy protons with the ra- 
diation field produced by the bright massive star in the 
system (e.g. a Be star with temperature T* ^ 3 x 10'* K 
in the case of LS I 303 and PSR B1259-63) occurs 

above Ep > [200 MeV/e,]mp ~ 2 x lO^^ [e,/10 eV] eV, 
where e, ~ 3fcr, is the typical photon energy of the stel- 
lar radiation. In contrast, high energy protons efhciently 
interact with the protons from the dense stellar wind al- 
ready at much lower energies. 

In the case of massive stars of type Be, the pp in- 
teraction rate can be highly enhanced if the high en- 
ergy protons are able to penetrate the dense equato- 
rial disk known to surround them. An obstacle for 
the penetration of the high energy protons accelerated 
e.g. close to the compact object into the disk could 



be the presence of magnetic fields, which would devi- 
ate the proton trajectories away from the disk. How- 
ever, the Larmor radius of the highest energy protons, 
Rl ~ 4x10^2 [Ep/10^^ eV] [B/1 G]"^ cm, could be com- 
parable to the size of the system. Thus, if the magnetic 
field in the region of contact between the stellar wind and 
the relativistic outflow is not larger than several Gauss, 
the highest energy protons can freely penetrate the dense 
stellar wind region. Note that the same effect may also 
result in the production of very hard injection spectra 
of the high energy protons in the pp interaction region 
(penetration of the lower energy protons into the dense 
stellar wind suppressed by the magnetic fields). 

The estimation of the pp interaction rate in the disk 
is straighforward, provided a known disk density pro- 
file. For example, in the case of the Be star in the 
LS I +61° 303 system, the density profile of the equatorial 
disc has been measured, using the IR free-free and free- 
bound radiation (see e.g. Refs. Waters et al. [2^, Martf 
& Paredes [1^). Assuming a disc half-opening angle of 
6*0 ~ 15° and a stellar radius of i?* = 10i?Oj the disk 
density profile can be modelled as a power law of the 
distance D 



nd(-D) ~ rid,o 



(2) 



where 7id,o ^ 10 
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cm~'^ and 7 ~ 3. 
Given Upp ~ 10"^^ cm^, the pp interaction cross- 
section at Ep ^ 10^^ eV, one finds the characteristic time 
of pp interactions in the disk 



tppiD) = 



1 



30 



CppTlc^ 



D 

rI 



1 T 



(3) 



To estimate the efficiency of the pp interactions, the pp 
interaction time should be compared to the escape time 
from the system. Since the details of the escape regime 
are not known, one can only compare ([3]) with possible 
relevant time scales. The shortest possible time scale is 
set up by the light-crossing time of the system. 
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20 



D 



(4) 



This time scale is relevant e.g. if the Larmor radius of 
the high energy protons is comparable to the size of the 
disk, so that these protons cannot be trapped in the disk 
interior. If the high energy protons escape on this short- 
est time scale, the pp interactions can only be efficient in 
the direct vicinity of the Be star, at D ~ 10^^ cm. 

The escape of the high energy particles can be slowed 
down in the presence of a sufficiently strong tangled mag- 
netic field B. In this case, the characteristic escape time 
is given by the diffusion time which is, in the Bohm dif- 
fusion approximation 
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FIG. 1: Mechanism of production of neutrinos in interactions 
of high energy protons ejected by the compact object with 
the dense equatorial disk of Be star. 



Even if the high energy particles are trapped by the tan- 
gled magnetic fields in the stellar wind, they do not stay 
in the innermost dense stellar wind region for a long time: 
the high energy particles escape together with the stellar 
wind of velocity v^, so that the escape time is limited 
from above by 



7 X 10-^ 



D 



10^ cm s^i 



(6) 



In the equatorial disk of a Be star, the initial wind veloc- 
ity close to the surface of the star is believed to be slow 
with ~ 1 — 20 km/s close to the surface of the star 
[Hi m, dOl ■ The polar component of the wind is faster 
and less dense than the equatorial one. Both equato- 
rial and polar winds accelerate with distance, reaching 
asymptotically the velocity Vca — (1-2) x 10'^ km s~^. 
Comparing ([6]) with ([3]) one finds that the pp interac- 
tions in the Be star disk can efficiently transfer the power 
contained in the high energy protons into neutrinos and 
7-rays. 



C. Broad- band emission from pp interactions 

pp interactions result in the production of pions, which 
subsequently decay onto 7-rays, neutrinos and elec- 
trons/positrons. The neutrinos freely escape from the 
source. The electrons/positrons suffer synchrotron and 
inverse Compton (IC) energy loss and, in this way re- 
lease their energy through 7-ray emission. 

If the Larmor radius of the highest energy protons is 
comparable to the size of the disk, most of the pions are 
produced by the protons propagating toward the com- 
panion star. In this case the neutrino and 7-ray emission 
from the pion decays is expected to be anisotropic, with 
most of the neutrino flux emitted toward the massive 
star, as it is shown in Fig. [1] Conversely, the synchrotron 
and IC emission from the e+e" pairs is, most probably, 
isotropized at energies at which the radiative cooling time 
of electrons becomes longer than the period of gyration 
in the magnetic field. Differing anisotropy patterns for 
neutrino emission and for broad-band emission from the 



secondary e+e" pairs should, in principle, lead to sig- 
nificant differences in the expected orbital lightcurves of 
neutrino and electromagnetic emission from the source. 

The IC cooling time of c+e~ pairs is shortest at the 
boundary between the Thomson and Klein-Nishina limits 
for Compton scattering. 3A:T*i? ^ (moC^)^, when ~ 30 
GeV: 



tic = 



37rm2c4£»2 



(T'tL-),E 

■1038 



2.5 



ergs" 



D 



1 2 



10^" eV 
E 



s, (7) 



At the same time, the period of rotation 
around the Larmor circle is = 2'kRl/c ~ 

1.3 [E/IQ^"^ eV] [B/l G]"^ ms< tic. 

The IC scattering on higher energy electrons proceeds 
in the Klein-Nishina regime. In this regime, e E, 
and the electron energy loss time grows with energy just 
slightly slower than the Larmor radius [3ll | 
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where the density of the blackbody photons has been 
diluted by the {D/Ri,)^ factor. This means that if the 
dominant mechanism of radiative cooling is the IC en- 
ergy loss, the emission from the secondary c'''e~ pairs is 
isotropic. 

The synchrotron cooling time, 
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can be much shorter than the period of rotation around 
the Larmor circle only at the highest energies 



E > 



{arBy/^ 



6 X 10^^^ 



B 



1 G 



-1/2 



eV 



(10) 



Taking into account that the synchrotron emission from 
electrons of energy E is emitted in the energy band 



eBE^ 
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1013 eV 



MeV, (11) 



this emission is expected to be isotropic below several 
MeV. 

Comparing the synchrotron and IC cooling times at 
energies £; < 30 GeV and iJ > 30 GeV, one finds that 
the lower energy electrons predominantly dissipate their 
energy via inverse Compton in the energy band 



eic ^ ikT^ 



E 



3 X 104 K 



E 



101" cV 



GeV, 
(12) 
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while synchrotron radiation losses dominate for the high- 
est energy electrons. 

A significant contribution to the cooling of electrons 
and positrons in the densest part of the pp interaction re- 
gion may arise from Bremsstrahlung emission, with cool- 



ing time ^Br 



10^ [7id/10^'^ cm] s, which compares 



to the synchrotron / IC cooling time of ^TeV electrons. 



D. Absorption of 10 GeV-TeV 7-rays 

The flux of 7-rays from the pp interaction region is 
absorbed due to the pair production on the UV photon 
background in the vicinity of the Be star. Maximal opti- 
cal depth with respect to the pair production is achieved 
at energies ~ AmUe^ ~ 0.2 [T*/3 x lO" K] TeV, 
where the pair production cross section reaches its max- 
imum a~f^ ~ 1.5 X 10^^^ cm^. The density of pho- 
tons close to the surface of the Be star is nph ~ 5 x 

l3 rr^ 1-2 



10" [r*/3 X 10^] [D/R^] cm-3, thus an estimate for 
the optical depth for 7-rays of this energy is given by 
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(13) 



assuming the radius of the star i?* ~ 5 x 10^^ cm. At 
higher 7-ray energies, E^T^, ^ (moC^)^, the pair pro- 
duction cross-section and the optical depth decrease as 
E'HnE. The attenuation of the 7-ray flux due to the 
pair production becomes small only at energies above 
~ 10 TeV. Below this energy, the spectrum of the 7- 
ray emission can be significantly different from that of 
the neutrino emission, so that no solid prediction for the 
neutrino flux can be derived based on the observed 7-ray 
flux and spectrum in the TeV band. 

The power of the absorbed 7-rays is re-distributed to 
the secondary e+e~ pairs, which subsequently loose their 
energy through synchrotron and inverse Compton emis- 
sion. Depending on the magnetic field strength in the 
pair production region, the bulk of the electromagnetic 
emission from the secondary pairs with energies ranging 
approximately between 10 GeV and 10 TeV can be ei- 
ther re-emitted back in the GeV- TeV energy band, if the 
inverse Compton loss dominates, or in the X-ray band, 
in case the synchrotron losses dominate (see Eqs. (fTTj) . 



III. NUMERICAL MODELLING OF THE 
BROAD-BAND SPECTRUM 



wind protons in the course of their escape from the sys- 
tem. The equatorial wind of the massive star is supposed 
to have the radial density profile described by Eq. ([2]). 
We assume that the high energy protons are initially in- 
jected at a distance Dmin — 1.2i?* and then escape to- 
gether with the secondary particles produced in pp inter- 
actions toward larger distances. We explore the different 
options for the escape regime, as described in Section [Til 

One possibility occurs when the high energy protons 
are traversing the pp interaction region without being 
trapped by the magnetic fields, the escape velocity is 
comparable to the speed of light. A similar "fast escape" 
situation is present when the magnetic field in the con- 
tact region between the stellar wind and the relativistic 
outflow is ordered (the situation present e.g. in the sce- 
nario of interaction of relativistic pulsar wind with the 
Be stellar wind). To model this "fast escape" regime, 
we assume that the primary high energy protons and the 
secondary e'''e~ pairs move toward larger distances with 
speed Wosc ~ c. 

Another possibility is that the high energy protons, 
which penetrate into the stellar wind, and/or the sec- 
ondary e+e" pairs, deposited throughout the stellar wind 
in result of pp interactions, are trapped by the tangled 
magnetic field in the stellar wind. In this case they would 
escape with velocity Vesc ~ Vw{D) approximately equal to 
the stellar wind velocity. The slow down of escape mostly 
affects the low energy parts of the electron spectra, at en- 
ergies at which the radiative cooling time is comparable 
or larger than the escape time from the system. 

We calculate the production spectra of neutrinos, 7- 
rays and e+e~ pairs using the approximations given in 
the Ref. 34]. The synchrotron, IC and Bremsstrahlung 
emission from the secondary e^e~ pairs, as well as the 
evolution of the spectra of the pairs due to the radiative 
cooling effects is modeled in a standard way [3l|, via a 
solution of the kinetic equations with the derivative over 
the distance Vcscd/dD substituted for the time derivative 
d/dt. 

To model the synchrotron emission, we assume a 
certain radial profile of the magnetic field, B = 
Bq (D/R^) with as = 1. The parameters Bq and 
as, used for numerical calculations shown in Figs. O [3] 
and 0] is Bq ~ 5 G, while for the calculation of Fig. [51 
Bq = 0.5 G. Taking into account that the trajectories of 
the secondary e+e~ pairs are isotropized by the magnetic 
fields and that we are interested in the orbit-averaged 
spectrum of the source, we take the angle-averaged cross- 
section of IC scattering for the calculation of the IC emis- 
sion. 



In order to demonstrate the possibility of dramatically 
different spectra of neutrino and 7-ray emission from the 
system in the TeV energy band, we have modeled the 
emission from pp interactions, assuming different proton 
injection spectra. Our numerical code follows the evo- 
lution of the secondary particle spectra produced by the 
interaction of the high energy protons with the stellar 



A. Broad-band spectrum of a 7-ray-loud binary in 
hadronic interaction model 

The broad-band spectra of the three known persistent 
or periodic 7-ray loud binaries, LSI 4-61 303, LS 5039 
and PSR B1259-63 have similar shape (which does not 
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FIG. 2: Broad-band spectrum of emission from secondary 
e^e~ pairs produced in pp interactions close to the surface of 
the Be star, calculated assuming a nearly monochromatic pro- 



ton injection spectrum with Tp — 0, i?p,( 



10'" eV. The up- 



per panel shows the injection spectrum of e e~ pairs (dashed 
line) and the spectrum formed as a result of cooling via syn- 
chrotron, IC and Bremsstrahlung emission as well as Coulomb 
energy loss. In the lower panel, thin red solid, dashed and 
dotted lines show respectively the synchrotron, the IC and 
the Bremsstrahlung emission from the pairs. The black thick 
solid line shows the overall broad-band model spectrum. 



resemble the shape of typical accretion-powered X-ray 
binaries [33|). The 7-ray-loud binaries spectral energy 
distribution are peaked in the MeV-GeV energy band. 
The observed X-ray and TeV emission apparently form 
the low and high energy tails of the MeV-GeV " bump" 
in the spectrum. Fig. [2] shows an example of LSI -1-61 
303 SED. 

Taking into account the possible anisotropy of the neu- 
trino and neutral pion decay emission, we first show in 
Fig. [2] the broad-band isotropic emission spectrum from 
the secondary e+e" pairs. The primary proton injection 
spectrum is assumed to be hard, with Tp ~ 0, so that 
most of the protons have the energy close to the cut- 
off energy assumed to be E'cut = 10^^ eV. As discussed 
above, such an almost monochromatic spectrum of pro- 
tons injected into the stellar wind can be produced when 
the high energy protons originate from a " cold" rclativis- 
tic wind with bulk Lorentz factor ~ 10^ or when only the 
highest energy protons have large enough Larmor radii 
to be able to penetrate deep into the stellar wind. 

Within the hadronic model of activity, the SED MeV- 
GeV bump can be ascribed to the synchrotron emission 
from the e+e" pairs produced in the pp and 77 inter- 
actions. One should note that if the primary proton 
injection spectrum is hard, the shape of the MeV-GeV 



synchrotron bump in the spectrum depends slightly on 
the details of the proton injection spectrum, because the 
shape of the e+c~ pair spectrum (shown in the upper 
panel of Fig. ^ is determined by the radiative cool- 
ing effects, rather than by the details of the e+e~ injec- 
tion spectrum (shown by the dashed line in the upper 
panel of the Figure). The hard e"^e~ pair spectrum be- 
low E ~ 10* eV is explained by the dominance of the 
Coulomb loss in the densest innermost part of the stel- 
lar wind. The E~'^ type spectrum between 10* eV and 
10^° eV is explained by the dominance of the synchrotron 
loss in this energy range. The hardening of the spectrum 
in the range 10^° eV < E < 10^^ eV is explained by the 
dominance of the IC energy loss proceeding in the Klein- 
Nishina regime. Above 10^^ eV, the IC cooling becomes 
less efficient than the synchrotron cooling, which leads to 
the softening of the spectrum. 

The electron energy at which the synchrotron loss takes 
over the Klein-Nishina energy loss can be estimated by 
comparing the synchrotron cooling time ^ to the inverse 
Compton cooling time ([8]): 



E 



KN-S 
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1 G 



cV 



(14) 



(where Bq is the magnetic field in the innermost part of 
the stellar wind and we assume the radial profile B ^ 
D^^). The synchrotron emission by electrons of such 
energies is produced in the energy band above 



€ ^ ^synch— bump — 3 
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McV, (15) 



resulting in a broad " bump" in the synchrotron spectrum 
above ckn-s- 

The comparison of synchrotron, inverse Compton 
and Bremsstrahlung emission spectra from the " cooling- 
shaped" electron spectrum to the LS I +61° 303 multi- 
wavelength data [IJ], done in the lower panel of Fig. [21 
shows that the model in which all the e"*"e~ pairs are 
initially injected at very high energies (above 100 TeV) 
provides a good fit to the broad-band spectrum of the 
source. To produce this figure, we have taken into ac- 
count the attenuation of the spectrum of the 7-ray emis- 
sion due to the 77 pair production in the photon field of 
the Be star, assuming that the 7-rays escape almost ra- 
dially from the emission region. This corresponds to the 
situation of the compact object close to the inferior con- 
junction of the orbit, where the effects of the absorption 
of the 7-ray flux arc minimized. 

Fig. [3] shows an opposite situation in which the 7-ray 
spectrum strongly absorbed in the energy band 0.1-1 TeV 
(long-dashed thick solid line) due to the pair production. 
This situation can occur when the emission comes from 
the densest innermost part of the stellar wind in the vicin- 
ity of the Be star. The real value of the average T~^.y can- 
not be estimated unless the details of the 3-dimcnsional 
geometry of the emission region and the relative orien- 
tations of the extended emission region, of the Be star 
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FIG. 3: Broad-band spectrum of emission from pp interac- 
tions, calculated assuming the same parameters as in Fig. 
[2l but considering the possibility of strong neutrino emission 
from the source. The notations for the 7-ray spectrum are 
the same as in the lower panel of Fig. |2l The green thick 
solid line shows the neutrino spectrum. Blue thick solid line 
shows the tt" decay contribution to the 7-ray spectrum. The 
black thick line shows the overall broad-band model spectrum. 
The dashed part of the line shows the part of the spectrum 
attenuated by 77 pair production. 



and of the observer are known. Taking into account this 
uncertainty, we choose for t~^^ the minimal value which 
ensures the non violation of the flux upper bound at ener- 
gies above ~ {severaT] TeV from VERITAS observations 
of the source [3J| . 

Fig. [3] also shows the 7-ray emission spectrum from 
neutral pion decays (blue thick solid line) and the neu- 
trino spectrum (green thick solid line). The observation 
of the MeV-GeV band synchrotron emission from the 
secondary e"'"e~ pairs enables to constrain of the possi- 
ble neutrino emission as explained previously. Contrary 
to the emission from e"'"e~ pairs, neutrino and -n^ de- 
cay emission is anisotropic. In principle, the anisotropy 
could lead to a boosting of the neutrino and tt*^ decay 
components. However, if the orbital plane is not aligned 
with the line of sight, the neutrinos and 7-rays from tfi 
decay should be emitted into a cone with rather wide 
opening angle, to be observable (the opening angle of 
the cone should be larger than the inclination angle of 
the binary orbit, see Fig. [1]). This means that the flux 
is at most moderately enhanced by the anisotropy effect. 
At the same time, if the inclination of the orbit is larger 
than the opening angle of the neutrino and tt*^ decay 7- 
ray emission cone, the flux in these emission components 
could be negligible. 

Figs, m and O show the results of calculation of the 
broad-band emission spectrum from pp interactions ob- 
tained assuming softer spectra for the primary protons, 
with Fp = 1 and Fp = 2, but keeping the same cut-off 
energy, E^nt = 10^^ eV. In these figures, we illustrate 
the effect of an additional uncertainty in the calculation 
of the broad-band spectrum of the source, related to the 
uncertainty of the cascade contribution to the observed 
source spectrum. The problem is that the shape of the 
electron spectrum below 10 TeV is modified by the in- 




log E [cV] 

FIG. 4: Same as in Fig. 3, but for Fp = 1. The thick short 
dashed line shows the overall spectrum calculated taking into 
account the tertiary e^e~ pairs from 77 interactions. 




log E [cV] 

FIG. 5: Same as in Fig. 4 but for the proton injection spec- 
trum Fp = 2. The magnetic field is assumed to be Bo = 0.5 G. 

jection of e^e~ pairs via the process of absorption of 
the 7-rays on the UV photon background. In case the 
synchrotron radiation dominates the energy loss of the 
"tertiary" electrons, the synchrotron emission from the 
cascade e+e~ pairs will further increases the height of the 
MeV-GeV synchrotron bump in the spectrum. In case 
inverse Compton dominates the energy loss of the cas- 
cade electrons, the tertiary pairs will emit mostly in the 
1-10 GeV energy band. The overall flux of the cascade 
contribution is determined by the average optical depth 
of the extended source with respect to the pair produc- 
tion. The thick short dashed lines in Figs. U] and O show 
the calculated modification of the broad-band spectrum 
of the source by the cascade. 

IV. ESTIMATE OF THE NUMBER OF 
NEUTRINO EVENTS FOR ICECUBE 

In the previous section we have shown that the 
hadronic model of activity of 7-ray-loud binaries could 
provide a good fit to the broad-band spectra of these 
sources. Within this model, the observed MeV-GeV 
bump in the spectral energy distribution is explained by 
the synchrotron emission from the secondary e'''e~ pairs 
produced either in pp interactions or in result of 77 pair 
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production process. 

A direct test of the hadronic model would be the de- 
tection of neutrino flux from the known 7-ray-loud bina- 
ries. In principle, the expected neutrino flux from a given 
source can be estimated from the known 7-ray flux. In 
the simplest case when the source is transparent to the 
7-rays, the estimate is straightforward, because the nor- 
malization and spectrum of the 7-ray flux produced by 
the neutral pion decays in the source directly gives an 
estimate of the spectrum and normalization of the neu- 
trino flux. However, modeling of the broad-band spec- 
tra of emission from 7-ray-loud binaries, described in the 
previous section, shows that this simple possibility, most 
probably, does not hold. First, the spectrum of 7-ray 
emission in the TeV band is affected by the 77 pair pro- 
duction, so that neither the overall flux, nor the spectral 
index of the 7-ray emission are related to the ones of the 
neutrino spectrum. Next, a sizable contribution to the 
spectrum of 7-ray emission in the TcV band can be given 
by the IC emission from the secondary e+e" pairs pro- 
duced in pp interactions. This further modifies the flux 
and the spectrum of TeV 7-ray emission. 

Thus, the information on the properties of the TeV 7- 
ray emission from a 7-ray-loud binary can not be used 
to estimate the neutrino flux from the source. Surpris- 
ingly, the source flux in the MeV-GeV, rather than TeV 
energy range can be used for the neutrino flux predic- 
tion. This possibility arises because within the hadronic 
model of activity, the MeV-GeV bump in the spectral 
energy distribution is produced by the emission from the 
secondary e"''e~ pairs from the pp interactions. Contrary 
to the TeV flux, the MeV-GeV flux from the 7-ray-loud 
binary is not affected by the pair production and can be 
used to estimate the total energy output from the pp in- 
teractions in the source. The only uncertainty of such 
an estimate is that the synchrotron emission from the 
secondary e+e~ pairs in the MeV-GeV energy band is, 
most probably, isotropic, while the neutrino and ir^ de- 
cay 7-ray emissions are not. As explained in the previous 
section, the observed source flux in the MeV-GeV band 
gives, in fact, an upper limit on the possible neutrino flux 
from the source. 

Although the total power of the neutrino emission can 
be estimated from the MeV-GeV luminosity of the 7-ray- 
loud binary, the modelling of the broad-band emission 
spectrum of the source gives only mild constraints on the 
neutrino emission spectrum: acceptable models of the 
broad-band spectra can be found assuming initial proton 
injection spectra ranging from a E~'^ power law to an 
almost monochromatic injection spectrum. 

Below we estimate the maximal neutrino event rate in 
ICECUBE for a particular example of LSI +61 303, tak- 
ing into account the uncertainty of the spectrum of neu- 
trino emission from the source. In the most optimistic 
case scenario, when the neutrino flux is at the level of 
the observed MeV-GeV flux from the source, we also in- 
vestigate the possibility to measure the parameters of the 
neutrino spectrum (and, therefore, of the primary proton 



spectrum), such as the spectral slope and/or the cut-off 
energy. 

The rate of neutrino-induced muon events from a point 
source at declination 6 with differential neutrino spec- 
trum d<i>^{E^)/dE^ 



= 



dE, 



J dE,Af{E,,S) 



(16) 
(17) 



where Acs{Ey^5) and A°'^(£^^) are respectively the neu- 
trino and the muon effective areas of the detector. The 
neutrino effective area, often provided in papers in con- 
trast to the muon effective area, folds the detector ef- 
ficiency and the neutrino propagation effects. Therefore 
Eq. p6p is convenient for a direct calculation of the num- 
ber of detected neutrino-induced muon events in a spec- 
ified neutrino energy range. 

However, the neutrino energy is not an experimental 
observable, contrary to the reconstructed muon energy. If 
we are interested in understanding whether the different 
neutrino spectra can be experimentally differentiated, i.e. 
we want to assess the spectral response of the detector 
for each model, we should use instead Eq. p?)) . This 
cannot be done without further derivation of the muon 
effective area and knowledge of the details of the detector, 
namely the angular resolution (of the order of 1° [i^ for 
ICECUBE in its final 80 string configuration (IC80) and 
precisely accounted for in the following) and the energy 
resolution of the order criogE « 0.3 [46l. |47|. 

The differential neutrino induced muon spectrum 
at the detector, for a given source at location (5, 
d<I>^(ii'^, (S)/di?^, is obtained after propagation of neutri- 
nos up to the interaction point and further propagation 
of the muons to the detector. This can be sketched as 
follows: 



ycc+E^) 



E, 



(18) 



Here -B^g {E^, ) is the muon energy at the interac- 
tion vertex, given by E^g{E^) = (1 - ycc{Ey))E^, 
where yccAEy) is the mean interaction inelasticity, 
taken from [5l| (calculations below will adopt a sim- 
plifying hypothesis, considering an average inelasticity). 
Pdct{Efi, Eu) is the probability density for a muon pro- 
duced with energy E^g{Ey) to reach the detector with 
energy i?^. for a given distance between the interaction 
vertex and the detector, i?^ is calculated from E^g{Ey) 
by integrating the muon energy loss rate dE^/dX = 
— {a + bEfj), where the a, h are the e nerg y independent 
standard rock coefficients from Ref. [5J| and X is the 
grammage traversed by the muon on its way to the de- 
tector, 



Pdct{E^, Ey) 



~dX{E^,E^,{E,))/dE^ 
Rf,{Efj,g{E^)) 
1 1 
\n{l + E,,{E.)/e) (E^ + e) 



(19) 
, (20) 
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where e = a/6. The function pAet{E^,E^) is set to zero 
outside of the interval € [0, £^^„(i?^)]. X{E^,E^j_^) 
is the grammage function for a muon of initial energy 
and final energy E^^ and R^{E^„) = X(0, E^„) is its 
range. Finally, the function which enters 

the final stage of calculation sketched in (fT8|) character- 
izes the energy resolution of the detector and is discussed 
below. 

The differential muon flux at the detector, which enters 
equ. pT|) can be calculated via the propagation of the 
muon flux from the interaction vertex as 



dEa 



J dE^PdctiEfj,,E^) 



X p,^(E,,6)p,^,iE,)^^^^. (21) 

where Pint (i?i/) = NACTcciE^)R^iE^g) is the neutrino in- 
teraction probability in the vicinity of the detector (po- 
tentially producing a muon within the reach of the detec- 
tor) and ptr{E^,S) = exp {—Na<j{E^)X{S)) is the earth 
transmission probability of a neutrino arriving from dec- 
lination 6 {i.e. after crossing a grammage X{S)). Na 
is the Avogadro number, ccc and a are respectively the 
charged current and the total muon neutrino cross sec- 
tions, taken from [45| . 

The separability feature of the kernel Pdct{E^i, E^,), to- 
gether with the implicit assumption of the absence of 
muon propagation fluctuation allows us to analytically 
extract Af by substituting Eq. ^ into (fTT]) : 



AfiE,) 



d ( 

V Pdot [E^, , E^ )pint ( El, ) 



dE^,{E,)/dE, 



(22) 



The neutrino effective area A°j^{E^) without absorp- 
tion is extracted from the neutrino effective area aver- 
aged over the whole northern hemisphere, {A^,^ {E^,5)) s 
from ^43^, using the relation 



Af{E,) 



{Af{E,,5))s 



/°^Ptr(-B.,e-7r/2)d(cos0) 



The transmission probability ptr {Eu , <5) which appears in 
the denominator (numerically computed using the Pre- 
liminary Earth Model [i^l) strongly decreases with en- 
ergy. Therefore, to estimate the neutrino effective area 
at the declination of LS I -1-61° 303, the effect of earth 
absorption is accounted for according to 



AZ (-E^, ^Lsi+61303) = Pti(£^i/, (^i 



LSI+61 



303 



(24) 

It should be mentioned that A'^ is approximately, at 
the modest energy considered here, the effective area at 
zero declination, assuming an isotropic detector response 
and analyses uniform down to the horizon (of course, 
analyses proceed non completely uniformly w.t.r. the 



declination to account for the difficult background rejec- 
tion near the horizon. Nevertheless, lately, as can be 
seen e.g. in [37| . there is a remarkable uniformity of the 
sensitivity, down to very close to the horizon). 

The effective areas are shown for IC22 (the ICECUBE 
22 string configuration) and IC80, in Fig. [6l In the lower 
panel, the ICECUBE muon effective area derived in (50l | 
shows a remarkable agreement above ~1 TeV, while the 
behavior at lower energy is rather different: our calcula- 
tion exhibits a soft exponential decrease below ~1 TeV 
contrasting to the sharp cutoff at ~100 GeV. Consider- 
ing the upper panel, we notice the fast decrease of the 
neutrino effective area at LSI-l-61 303 declination above 
~300 TeV. This is due to the fact that the neutrino tra- 
jectories towards ICECUBE are slightly passing through 
the high density earth core. This leads to a signifi- 
cant suppression of the neutrino flux above this energy 
(the observation of PSR B1259-63 and LS 5039 with the 
Km3Net[3H| would not suffer this limitation and thus of- 
fer the possibility to probe spectral features up to the 
highest cutoff energies). 

Using the analytical muon effective area (P^ . the dif- 
fuse atmospheric neutrino background rate and the point 
source signal event rate as a function of the chosen recon- 
structed muon threshold energy i?thr can be calculated. 



7^b(^7 >Eti,,,^,5) = j dE, 



d(f>-'-(i?^,5) 



dE„ 



( 

xh{^)A''^{E^)u,,{E^,E,^,) (25) 



f d$™ 

7^,(£;;^^ > e,^,,^,5) = y di?^ — 



dE,. 



X E^) Af{E^)LO,, (S^, i?thr) 



(26) 



/b(V') ~ """0^ is the solid angle subtended within a radius 
Tp for the diffuse atmospheric neutrino background and 
fs{fp, E,j) is the fraction of signal events with energy Ey 
reconstructed within an angle ij) from the source: 



/s(V', E^)= ^ fli 1 - exp 



2^2 , 



, (27) 



(23) where the point spread function (PSF) cri,<p(£',y) = 
'* + ^liiiEu) decomposes into two distinct PSF 



contributions: the energy-dependent neutrino to muon 
kinematics Qi,^{E^) = 0.54°/^i?,y/TeV and the muon 
reconstruction resolution cri.*^,^,^^^ . For our analysis we 
extracted the parameters ai^qi and from a fit to the 
IC22 and IC80 PSF found in 1^ (for IC80. we obtain: 
flicso ^ 0.44, (jicso = 0.47°, 4^80 = o.51, tjicso = 1.12°. 
E.g. at 10 TeV, 50% of the neutrinos coming from a 
specified source are reconstructed within 0.9° for IC80). 
In the following, the results will be presented using the 
optimized = 1.3° for IC80. 

Above, in eas[25land[26l the energy argument for fs is 
set to E^, = E^ . This is a conservative statement because 
the muon energy is smaller than the neutrino energy, we 
have/,(V',i;j < U{-P,E,{E^)). 
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FIG. 6: Top: Comparison of different neutrino effective ar- 
eas. Red dashed line showns the effective area Af with no 
absoption. Blue dot-dashed line shows the effective area after 
averaging over the declinations 5 for an isotropic neutrino flux 
{Af) O. Curves running in parallel are for IC22 and IC80. 
Black solid line shows the effective area for LS I -|-61° 303. 
Bottom: Full lines are the muon effective area calculated 
with for IC22 and IC80. The dashed curve is the muon 
effective area derived in [soll. 



The weight factor used to calculate the number of 
events above some reconstructed energy: 



p,,(^;-^i?^)d^;-^ (28) 



follows from an assumed gaussian energy resolution 
distribution function in Alogi^, , ) = 

Va^ (log E'^^ - log E^), with GE = 0.3. 

In our calculations, the muon propagation was imple- 
mented through —{dE/dX) = a + bE. with the energy- 
independent parameters a and h for standard rock (ad- 
equate, given the source location and the proximity of 
ICECUBE to the bedrock), in order to derive the muon 
energy distribution at the detector [53|. Therefore, fluc- 
tuations of muon energy due to the stochastic nature of 
high energy muon energy losses which implies muon en- 
ergy distributions, were not accounted for (49| . A more 



refined treatment, which requires numerical methods, 
proceed in integrating over the distance the probability 
F{Ef^g, Efj^,l) [il, [HO] for a muon generated at distance I 
with energy E^^^ to reach the detector with energy £^^. 
Our simplified treatment has nevertheless the advantage 
of transparently combining the latest published neutrino 
effective area with some phenomcnological assumptions 
for the detector muon effective area, detector response, 
muon propagation and neutrino interaction, resulting in 
a simple analytical formula for the rates, thus avoiding 
more generic integro-differential forms. 

We choose the angular-dependent parametrization 
from Volkova (prompt -I- conventional neutrino fluxes) for 
the atmospheric neutrino background model [4l|. This 
corresponds to 7.7 detected muons when using the neu- 
trino effective area and per year with IC80 in the 1.3° 
circular bin at LS I -f61° 303 declination. Conversely, cal- 
culating this number with the analytical muon effective 
area [22] yields 8.8 detected muons. This 15% difference 
may come from the various simplifications and is proba- 
bly dominated by the delta function differential cross sec- 
tion approximation (at low energies, the average range of 
the neutrino induced muons is larger). At LS I -1-61° 303 
declinations, the atmospheric neutrino zenith angle to be 
considered is 9 ^ 5 - arcsin (2 sin^ 5 - 1) = 29° (6* = 0° 
is for vertical incidence). 

The results, in terms of the atmospheric background- 
subtracted muon spectra for the baseline neutrino flux 
models discussed in the previous sections, are presented 
in Fig. [71 The exposure time is taken to be 3 years of run- 
ning the full ICECUBE array. For comparison, we also 
show by the red solid thick line the 5a level above the 
atmospheric background (extracted from binomial statis- 
tics, where the trial factor is the total expected number 
of atmospheric neutrino n{E) in a 2'0 width declination 
band reconstructed with energy > E and the probabib- 



lity of success is p = 
tion band, nr^^{E) 



7r?/'^/27r Jg^f dcos0 in the declina- 



np + 5-\/ np{l — p). Note however 
that an unbinned calculation slightly lower the signal re- 
quirements fssji. 

We calculate the number of muon events as a function 
of the energy threshold for the three neutrino spectra 
discussed in the previous section. The blue dotted line 
show the detected muon spectrum for the model with the 
proton injection spectrum with spectral index Fp = 2, 
the solid black line corresponds to the proton injection 
spectrum with Fp = 1 and cut-off at the same energy, 
while the magenta dashed line corresponds to Fp — 0. 
In all three cases the cut-off energy is assumed to be 
-Bcut = 1 PeV. 

One can see that the uncertainty of the proton injection 
spectra results in a factor of ~ 2 uncertainty in predic- 
tions of the number of neutrinos detectable in ICECUBE. 
Assuming that the neutrino flux saturates the upper limit 
imposed by the observed 7-ray flux in the MeV-GeV en- 
ergy band, one flnds that the neutrino source should be 
detectable at 5cr level in the energy band 1 — 10 TeV in 
roughly one year, if the proton injection spectrum is not 
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too hard. 

Inspecting the muon spectra shown in the upper panel 
of Fig. [71 one sees that softed proton injection spectrum 
results in a slight excess of muon events at lower ener- 
gies. If the overall normalization of the neutrino flux 
were known, the measurement of the muon event spec- 
trum would allow to constrain the spectrum of the pri- 
mary protons in the source. However, taking into account 
the uncertainty of the overall normalization of neutrino 
flux introduced by the uncertainty of the anisotropy pat- 
tern of neutrino emission, one finds that the statistics of 
the signal will in the best case only allow for marginal 
estimate of the neutrino spectrum parameters. This is 
illustrated in the lower panel of Fig. [7l where a compar- 
ison of the shapes of the muon spectra is shown. If one 
assumes the same total number of muon events, the dif- 
ference in the spectra for the three models mostly remain 
within ^ la error bars over the entire energy range. 

The neutrino effective area in the direction of 
LS I +61° 303 peaks at the energy - 100 TeV (see Fig. 
ini Taking into account the typical inelasticity of pp in- 
teractions, one finds that neutrinos of this energy are in- 
jected by the primary protons with energies E'p ~ 1 PeV. 
This means that assuming the cut-off in the proton spec- 
trum at the energy -Ecut ^ 1 PeV, one maximizes the 
number of neutrino events detectable with ICECUBE. If 
the cut-off in the proton spectrum is much above 1 PeV, 
the neutrino signal will be lost because of the opacity of 
the Earth. Otherwise, if the cut-off in the proton spec- 
trum is much below 1 PeV, the signal is lost because of 
the transparency of the detector for the neutrinos and 
because of the high atmospheric neutrino background. 
Thus, the typical numbers of neutrino-induced events in 
the ICECUBE, shown in Fig. [7] should be considered 
as the upper bound on the possible neutrino signal from 
LS I +61° 303. 



V. CONCLUSIONS 

We have estimated the neutrino fiux from GRLBs ex- 
pected within the hadronic model of activity of these 
sources. Within such a model, the measured spectral 
characteristics of 7-ray emission from the source in the 
TeV energy band are not directly related to the spectral 
characteristics of the neutrino emission, because of ab- 
sorption of the TeV 7-rays on the thermal photon back- 
ground produced by the massive star in the system. The 
uncertainty of the calculation of the attenuation of the 
TeV 7-ray flux introduces a large uncertainty to the es- 
timate of the neutrino flux based on the measured TeV 
7-ray flux. 

Taking this uncertainty into account, we have adopted 
a different approach for the estimate of the neutrino flux 
from a GRLB. Namely, we have noted that the energy 
output of proton-proton interactions, and hence the neu- 
trino flux, can be constrained by the broad-band spec- 
trum of the source. The idea is that the e+e" pairs, pro- 
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FIG. 7: Background-subtracted cumulative muon spectra 
A^(i^/j,thr), expected after the 3-year ICECUBE exposure for 
the three model neutrino spectra of LS I +61° 303, discussed 
in the previous sections (error bars of signal muon spectra 
are the sum in quadrature of statistical errors of signal -|- 
atmospheric neutrino background). The blue dotted / black 
solid / magenta dashed line respectively show the spectra cor- 
responding to proton injection spectra with spectral indices 
r = 2, r = 1 and F = 0. The red thick solid line shows the 
5a excess above the atmospheric neutrino background (the 
" discovery threshold" ) . 



duced in the decays of charged pions, release their energy 
in the form of the synchrotron, IC and Bremsstrahlung 
emission, which contributes to the observed broad-band 
SED of the source. The requirement that the power of 
electromagnetic emission from the secondary c+c~ pairs 
in different energy bands is not higher than the observed 
source luminosity imposes a constraint on the neutrino 
luminosity of the GRLB. 

We have worked out the numerical model of the broad- 
band emission from a GRLB, assuming that pp interac- 
tions take place in the vicinity of the bright massive star, 
in the innermost and densest part of the stellar wind. 
We have shown that the interactions of high energy pro- 
tons (with energies reaching '--^PeV) in the dense stel- 
lar wind can explain the observed broad-band emission, 
with the MeV-GeV "bump" of the SED being due to the 
synchrotron emission from the secondary e"'"e~ pairs pro- 
duced either in the decay of charged pions or resulting 
from the absorption of 10 GeV-TeV 7-rays on the thermal 
photons from the massive star (see Figs. 1-4). 

Although the observed bolometric luminosity of the 
source (i.e. the height of the MeV-GeV bump of the SED) 
constrains the overall neutrino luminosity, the shape of 
the neutrino spectrum and the overall normalization of 
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the neutrino flux arc only mildly constrained by the 
multi-wavelength data. The problem is that the proper- 
ties of the broad-band emission from the secondary e+e" 
pairs are mostly determined by the effects of the radia- 
tive cooling and the escape of the pairs from the source, 
rather than by the initial injection spectrum of the pairs. 
This means that quite different injection spectra of the 
primary high energy protons (and hence of the neutrinos) 
can result in approximately the same broad-band SEDs, 
as is clear from Figs. 1-4. In addition, the anisotropy 
pattern of neutrino emission is expected to differ from 
that of the synchrotron emission from the e+e^ pairs. 
This results in an uncertainty of the anisotropy-related 
suppression of the neutrino flux compared to the MeV- 
GeV synchrotron flux. 

Taking into account these uncertainties of the neutrino 
emission spectrum, we have estimated the expected num- 
ber of neutrinos which will be detected by ICECUBE, as- 
suming that the neutrino flux saturates the upper bound 
imposed by the observed 7-ray flux in the MeV-GeV 
energy band. Considering the particular example of 
LS I -1-61° 303, we have found that if the spectrum of 
high energy protons in the source extends to the PcV 
energies, the source would be readily detectable within 



roughly one year of exposure with ICECUBE. 

We have also explored the potential of the full ICE- 
CUBE detector for the measurement of the spectral char- 
acteristics of the neutrino signal from LS I +61° 303. We 
find that in the case when the neutrino flux is at the level 
of the upper bound imposed by the observed MeV-GeV 
7-ray flux, an exposure time longer than 3 years will be 
required to reliably constrain the spectral index of the 
primary high energy proton spectrum via observations of 
neutrino signal in ICECUBE. The mere fact of a possible 
spectral characterization of neutrino sources within the 
lifetime of ICECUBE represents however a most exciting 
prospect for the future of mutimessenger astronomy with 
the possible emergence of a new field, the astrophysics of 
neutrinos. 
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